Topological Theory of Classical and Quantum Phase Transition 
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We presented the topological current of Ehrenfest definition of phase transition. It is shown that 
different topology of the configuration space corresponds to different phase transition, it is marked 
by the Euler number of the interaction potential. The two phases separated by the coexistence curve 
is assigned with different winding numbers of opposite sign. We also found an universal equation of 
coexistence curve, from which one can arrive the phase diagram of any order classical and quantum 
phase transition. The topological quantum phase transition theory is established, and is applied to 
the Bose-Hubbard model, the phase diagram of the first order quantum PT is in agreement with 
recent progress. 
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A. Introduction 

A phase transition(PT) brings about a sudden change 
of the macroscopic properties of a system while smoothly 
varying a physical parameter. While in topology, this 
sudden change is always marked by Euler number. For 
example, when the radii R of the mouth of a rubber bag 
smoothly shrinks to nothing, the Euler number jumps 
from 1( Euler number of a rubber bag) to 2(Euler number 
of a rubber sphere) at R — 0. One can not help asking: 
does PT bear a topological origin? 

On the origin of PT, Caiani et al proposed the topolog- 
ical hypothesis^] the PT is determined by a change 
in the topology of the configuration space, and the loss 
of analyticity in the thermodynamic observables is noth- 
ing but a consequence of such topological change. This 
hypothesis has received direct evident support from nu- 
merical and analytical computation of the Euler char- 
acteristic in 4 lattice modelQ, mean-field XY modelQ 
and k-trigonometric modelj^]. It is also supported by the 
Bishop-Peyard model of DNA thermal denaturationQ. 
Recently, it is proposed that the occurrence of a PT is 
connected to certain properties of the potential energy 
vQ. Franzosi, Pettini and Spinelli provided a rigorous 
theoretical support to this topological hypothesis 8] 9] . 

In this letter, we shall start from a more practical point 
to study this topological hypothesis of PT. It will be 
shown that the topological number of the critical points 
of the microscopic interaction potential naturally ap- 
peared in the topological current of Euler number. Fur- 
ther more, based on Ehrenfest definition of PT, we found 
an universal equation of coexistence curve, from which 
one can arrive the phase diagram of any order PT. This 



equation holds both for classical and quantum PT. From 
our definition of quantum PT, we derived the phase di- 
agram of superfluid-Mott-insulating quantum PT in the 
Bose-Hubbard model, it is in perfect agrement with the 
well known phase diagram. 



B. Topological theory of classical phase transition 

For a standard Hamiltonian system, which is described 
by H = XJipi \v\ + V(qi, ?jv), the Helmoltz free en- 
ergy readB0>G9) = -(2/3)- 1 L « ? ( 7 r//3) - /(/3)//3, where 
(3 = 1/T is the inverse of temperature. The first part 
of F{[3) comes from the kinetic energy and the second 
term /(/?) = jjLog J d N q exp[— f3V(q)) comes from the 
microscopic interaction potential V(q). 

As we know, a PT happens when we change the tem- 
perature or pressure of the system. In a more practical 
approach, we introduce two extra microscopic practical 
parameter {71,72}, which represent pressure, external 
magnetic field, electric field, and so on. The interac- 
tion potential {71,72} is generalized to ^(5,71,72)- In 
the frame work of Ehrenfest definition, the free energy 
F(ji, 72) is utterly depend on the two extra parameters. 

According to Ehrenfest definitional of a PT, when 
the system transforms from phase A to phase B, the 
first order phase transition requires that the free en- 
ergy ^(71,72) is continuous at the transition, i.e., 5F = 
F B (71,72) -F A (7i,72) = 0, but 



d 7l SF^0, d l2 SF^0 



(1) 



where d 7l — d l2 = For the second order 
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transition, it requires 



5F = 0, d 11 SF = 0, djJF = 



parameters {71, 72}, the Helmoltz free energy reads 
Fhi, 72) = -VP)" 1 Log(n/ (3)- /(/?, 7l , 72 )//3. with 
/(A 71, 72) = j f LogJd N qexp[-pV(q, ll , 72)]. Then 



d^JF 7^ 0, d^d^SF ^ 0, d 2 2 <5F ^ 0. (2) t h e gradient field of free energy is 



The pth-order PT is characterized by a discontinuity in 
the pth derivative of the free energy. Ehrenfest definition 
suggests that the change of the free energy between the 
two phases SF — F B (-fi, 72) — F A ('ji, 72) plays an essen- 
tial role in defining the order of PT. So we choose the SF 
as the fundamental order parameter field. 

The order parameter of the pth-order PT, 0, can be 
chosen as the (p — l)th derivative of SF, i.e., 1 — 
d 1 £^r l 8F 1 2 — d^T 1 SF. In order to investigate the topo- 
logical properties of a PT, we map the order parameter 
into a two dimensional sphere by introducing an unit 
vector field n with n a = a /||0||, n a n a = 1, where 
1 10| I = yt(f) a (f> a (a = 1,2). It can be proved that the 
Gaussian curvature on the two dimensional extra-space 
can be expressed in terms of ft: 



J = £ 

i ,j, a, b 
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In light of Duan's 
theory [ll|, using 9^ = 
Green function relation in 

&{Z\ (X _ d 



mapping topological current 

~ ^ ^ adi w\ and the 
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-space : d a d a In 
2irS 2 (4>), {d a — af^), we can prove that 
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£ 2 (0>(^)=5 2 (0){0\0 2 }, (4) 
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where D(cj>/q) = ^e jk e a bdj(f> a dk4> b is the Jacobian vec- 
tor, in the extra-two dimensional space, it is actually the 
Poisson bracket of 1 and 2 , {0\0 2 } = £<(§£££ " 

80 1 dtp 2 \ 

In Ehrenfest definition, the (p — l)th order PT re- 
quires that 4> a 7^ 0, (a = 1, 2), the pth order PT requires 
a = 0, (a = 1,2). Eq. (@J suggests that the topological 
current of Gaussian curvature: J = 0, iff 7^ 0; J ^ 
0, iff cf> = 0. So this topological current vanishes for 
the (p — l)th PT, and only exist for the pth order PT. 
In other words, the Euler number of the (p — l)th PT is 
zero, and it is not zero for the pth order PT. 

Therefore the topological information of the system 
focused on the zero points of the vector field 0. In fact, 
zero points of the vector field just corresponds to the 
extremal point of the variation of free energy, since the 
vector field is defined by the gradient field of SF. 

For a classical system, when the microscopic inter- 
action potential V(q) depends on the two microscopic 



1 Jd N qexp[-(3V}§^ 



OF 



<9 7l N J d N q exp[-(3V] 
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(5) 



Now we see the topological current relies on the criti- 
cal points of V in configuration space. Therefore it is 
the extremal point of the interaction potential that de- 
cides a PT. Eq. ||SJ also suggests that the higher order 
is decided by powers of a few basic ingredient, such as 
dV/d^/i, d 2 V/djid^fj, d 3 V/d^idjjd~fk, and so on. 

The implicit function theory shows, under the regular 
condition|lj F>(0/ 7 ) + 0, i.e., {0\ 2 } ^ 0, we can solve 
the equations = and derive n isolated solutions: Zk — 
(7fc)7feL (k = 1,2,... n). According to the S— function 
theory [l^J, one can expand S(4>) at these solutions, the 
topological current of Gaussian curvature is rewritten as 



J = J2^ r lkS(j-j 1 k )S( 1 - 1 2 ), 



(6) 



k=l 



here l3k = | | is the Hopf index, 
Wk is the winding number around Zk- 
rjk =signD(<f)/q) Zk =signdetM S F(zk) =sign{0\0 2 } 
■ 1 is the Brouwer degree, in which Msf{z%) is Hessian 
matrix. 

The winding number flkVk = Wk is the generalization 
of the Morse index in Morse theory. Its absolute value (3k 
measures the strength of the phase, rjk =sign{0 1 ,0 2 } = 
±1, represents different phases. For example, in the 
phase diagram of ultra cold atoms in an optical lat- 
tice(FIG. 1), in the area circulated by the curve, 
r]k =sign{0 1 ,0 2 } = +1, it represents the Mott- Insulator 
phase, while outside the curve, rjk =sign{0 1 , 2 } = —1, it 
represents the superfluid phase. On the curve, {0 1 , 2 } = 
0, that is where a quantum PT takes place. Loosely 
speaking, the PT which happened at the critical point 
with higher winding number is obviously different from 
the process happened at the point with lower winding 
number. 

For a more sophisticate case, when there are many 
phases in the phase diagram, the total winding number 
of different phases, X — J Jd 2 q — J2k=i W*> is the Euler 
number. There would be a substantial change of topol- 
ogy of the configuration space, when the phase diagram 
jumps to a higher order PT or a lower order PT. Such 
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as, the manifold can be jumps from a sphere to a torus, 
and to a torus with many holes, the Euler number would 
jumps from \ — 2 to \ — 0> then to \ = 2(1 — h) with 
h as the number of holes of the torus. During these pro- 
cess, the phase transition jumps from first order to the 
second order, and to the pth order. 

C. The universal equation of coexistence curve 

The most important thing we concerns about is the 
phase diagram of a PT. If we have found the coexistence 
curve equation, the phase diagram can be derived by solv- 
ing this equation. Our topological theory provide us an 
universal coexistence curve equation: 



{0\0 2 } = O. 



(7) 



It is an unification of the special coexistence equations of 
different order PT. As all know, in quantum mechanics, 
if two operator A and B commutate with each other, i.e., 
[A, B] — 0, they share the same eigenfunction. So Eq. Q 
reflects the same properties of different phase operators 
cj) 1 and cj) 2 . 

In the bifurcation theory of topological current 0], the 
topological current will split into several branches, when 
D(cj)/q) = 0, i.e., <j> 2 } = 0. It is rather surprising that 
the coexistence curve in the phase diagram happened to 
be the bifurcation curve of topological current. 

We first take the second order PT for example to 
demonstrate this. The order parameter of the second 
order PT is (j) 1 — drSF and cj) 2 = dpSF, substituting 
them into the Jacobian vector 



= ww- W W = ' 

and using the relations 



(8) 



OtOtSF = 



T 



dpdpSF = V(k$ - k|), 



d P d T 5F = V(a B -a A ), 



(9) 



we arrive 



V 



DWq) = -(C»-C*)( K %-^)-(Va"-Va A ). (10) 
Recalling the Ehrenfest equations 

dP _ a B - a A dP _ C B - C A 

dT k b -k a ' dT TV(a B - a A ) ' 1 ' 

it is easy to verify that the equation above is in consistent 
with the bifurcation condition Eq. (|l(Jf> . So the bifurca- 
tion equation D(<j>/q) = is an equivalent expression of 
the coexistence curve equation. 



For the first order PT, we chose the vector order pa- 
rameter as 4> — d°SF, here '0' means no derivative of the 
free energy. The generalized Jacobian vector of the first 
order PT with <fi = d°5F is given by 

^,6, ,dF B dF A , f dF B dF A n , . 



in mind of the relation 



i)F 

or 



-S and §f 



considering D(&) = 0, we have 



dP 
dT 



{S B - S A ) 
(V B - V A ) ' 



V, and 



(13) 



This is just the famous Clapeyron equation. 

The coexistence equation {(f) 1 , <fi 2 } = also holds for 
a higher-order PT. We consider a system whose free 
energy is a function of temperature T and magnetic 
field B, then the Clausius-Clapeyron equation becomes 
dB/dT = — AS/ AM. If the entropy and the magne- 
tization are continuous across the phase boundary, the 
transition is of higher order. For the pth order PT, the 
vector field is chosen as the (p — l)th derivative of SF, 



cf, 1 = gP^SF, (j) 2 = d^SF. Substituting 
Eq. @, we arrive 



D{<t>/q) 



d p 5Fd p 5F ddP^SF dd^SF 
OTP dBP ~ 8B8TP- 1 dTdBP- 1 

d 2 F 



into 



0. (14) 



Considering the heat capacity 



-Hp- and the sus- 
ceptibility t^!- = x, the bifurcation condition D{<p/q) = 
is rewritten as 

~dB 



dT 



(-iy 



Ad p ~ 2 C/dT 



•p-2 



(15) 



T c AdP~ 2 x/dBP- 2 ' 

This is in excellent agreement with Kunmar's results 

So different phases coexists on the the bifurcation 
curve, when PT occurs, the system bifurcate into several 
phases, each of them is assigned with a winding num- 
ber, but the total topological number is conserved, it is 
a constant which relies on the topology of configuration 
space. 

In fact, the Jacobian vector also provide us a new spec- 
tacles to the scaling lawfl^. As all known, at the crit- 
ical temperature, the correlation length £ becomes infi- 
nite, from which the divergence of the physical param- 
eters arises The topological current is invariant 
under renormalization group transformation q' = R(q). 
The new topological current J' expressed by the new pa- 
rameter q' = (q 1 ,q 2 ') share the same fixed point with 
J, J' = J = oo. This is an open question which need 
further study. 
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D. Topological quantum phase transition theory 

So far there is still no such a general topological theory 
of quantum PT, in the frame work of this paper, it is easy 
to establish an unified topological quantum PT theory. 

The usual thermodynamic phase transition at finite 
temperature is driven by thermal fluctuations, as tem- 
perature is lowered, the thermal fluctuations are sup- 
pressed. But the quantum fluctuations still exist and 
play a vital role in driving the transition from one phase 
to another. At the absolute zero of temperature 0], 
quantum PT can be accomplished by changing not the 
temperature, but some parameter in the Hamiltonian of 
the system. This parameter might be the charging en- 
ergy in Josephson-junction arrays (which controls their 
superconductor-insulator transition), the magnetic field 
in a quantum-Hall sample, doping in the parent com- 
pound of a high-Tc superconductor, or disorder in a con- 
ductor near its metal-insulator transition. 

The order parameter field of a quantum PT can be 
chosen as the ground state energy! 



E g = —In Z, Z = J D<j>e 

It plays a similar role as free energy 
1 



i J dxdtC(4>) 



F 



-In Z, 



Tr e - /3ff(7i( * )) . 



(16) 



(17) 



In fact, they are equivalent, we considering a hermitian 
Hamiltonian H(t) — H(~fi(t)) which is time-dependent 
through a set of parameters 7i(t), i = 1, n. Suppose, 
that for any fixed t on the time interval [0, T], the spec- 
trum of H is non-degenerate: if(7)|n(7)) = £' n (7)|n(7)), 
with 71(7)) as an eigenstate of Hamiltonian. 

As seen above, the gradient field of the free energy 
plays the keys role in a classical PT, here we also choose 
the order parameter of quantum PT as 

dF i Tre-^W)gf 
<9 7i /3 Tr e -/3H( 7i (t)) 
In mind of Feynman-Hellman theorem: 



(IS) 



one can sees that the 



(n(7)I^K7)> = *f£l, 
order parameter of quantum PT is actually the gradient 
field of the variation of the ground state energy 8E g . 

In analogy with Ehrenfest definition, the first or- 
der quantum PT, in which the system transforms from 
phase A to phase B, requires that the ground state en- 
ergy E g (71, 72) is continuous at the transition, SE g = 
Ff( 7l ,7 2 )- ^(71,72) =0, and 



where <9 7l = d l2 
quantum PT, it requires 



For the second order 

072 



SE g = 0, d^SEg = 0, d l2 5E g = 0, 
d^SE g ^0, d 71 d l2 5E g 0, d^SEg^O. (20) 

The pth-order quantum PT is characterized by a discon- 
tinuity in the pth derivative of difference of ground state 
energy SE g , 

d^SEg = 0, d^d^-HEg = 0, (m = 1, 2, ...,p- 1), 
d^HEg = 0. 

d^SEg + 0, d™cF- m 5E g + 0, (m = 1, 2, ...,p), 
d p 72 SE g £ 0. (21) 

In the topological quantum PT theory, we can chose 
the vector order parameter (\> as the (p — l)th derivative 
of SE„, 



1 = dP-HE g , 2 = dP- 1 5E : 



72 



(22) 



Then the phase diagram can be obtained from the coex- 
istence curve equation, 



f.1 PlJ,2 



971 dj2 ^72 371 



= 0. 



(23) 



d^SEg ytO, d^SEg^O, 



(19) 



Here 7^ is the physical parameters of a quantum system. 

A gas of ultra cold atoms in an optical lattice 
has provided us a very good experimental observation 
of superfluid-Mott-insulator phase transition 20] . In 
gapped quantum spin systems, an explicit divergence of 
the entanglement length appears in the ground states j^] 
at a topological phase transition. Using Matrix Prod- 
ucts States, it was found that the localizable entangle- 
ment possess a discontinuity at its first derivative at a 
Kosterlitz-Thouless phase transition, 22]. In the follow- 
ing, we shall establish a complete theory of topological 
quantum PT. 

This system is described by the Bose-Hubbard model 

H = - J ^2(a]aj +a]a,j) +^£i7"ii + y^,ni(nj - 1), 

(24) 

the last term U corresponding to the on site repulsion 
between atoms, while the first term J describes the tun- 
nelling of atoms between neighboring sites. 

Using mean-field approaches, the ground state energy 
of Bose-Hubbard model can be generally expressed as 
the functional of £j, J and U, i.e., E g = E g (si, J,U). 
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here we have taken 71 = U, 72 = J- From this equation, 
one can arrive the phase diagram of the pth order quan- 
tum PT from the equation above, and find some new 
quantum phases. 



FIG. 1: The phase diagram of the Bose-Hubbard Hamiltonian 
for g=l. Inside the curve is the Mott-insulating phase, outside 
the curve is the superfluid phase. 



Following perturbation theory up to the second order, 
the variation of ground state energyji^ SE g is 



SE g = [ 



1 



U(g - 1) - J J — Ug 



(25) 



From our definition of the first order phase transition, the 
boundary between the superfluid and the Mott insulator 
phases should be decided by the equation 5E g = 0. We 
have plotted the phase diagram(FIG. 1). This phase dia- 
gram is in perfect agrement with the well-known phase di- 
agram of superfluid- Mott-insulator phase transition)^. 
It is also the same as the phase boundary obtained from 
the minimizing the free energy [2^, when treat the hop- 
ping term as perturbation. So most of the recent phase 
diagram of the quantum phase transition are focused on 
the first order. 

As for the pth-order quantum PT, it is characterized 
by a discontinuity in the pth derivative of the variation 
of ground state energy. The order parameter field of 
the pth order Quantum PT can be chosen as the (p- 
l)th derivative of 6E g with respect to U and J, i.e., 
^ = cfj- x 6E g , f^^SEg. 

The equation of coexistence curve is given by 
{4>\^} =D{±) = 0, 



xv v f dJP dUP 



ddP 



^5 E g ddP~ 



l 5E n 



dUd.JP- 1 dJdUP- 



E. conclusion 

In summary, when a PT jumps from lower order to 
higher order, the Euler number of the extra configura- 
tion space would jump from one integer to another. In 
other word, different order of PT corresponds to different 
manifold with different Euler number. 

We have found an universal equation of coexistence 
curve equation, {(j) 1 , <p 2 } = 0, from which one can derive 
the phase diagram of any order. In the phase diagram, 
different phase are assigned with different winding num- 
bers, the sum of these winding numbers is Euler number 
which bifurcated on the coexistence curve. 

In analogy with Ehrenfest definition of classical PT, 
we proposed the higher order quantum PT, in which the 
ground state energy correction 5E g plays the key role. 
The universal equation of coexistence curve equation, 
{(f) 1 , cj) 2 } = 0, also holds for quantum PT. 

Now we proposed the general process to a obtain the 
phase diagram of pth order quantum PT. First, we calcu- 
late the ground state energy correction using perturba- 
tion theory. Then, substitute SE g into {(j) 1 , <fi 2 } = 0, and 
solve this equation. This method is an universal theory, 
it can be generalized to a great variety of classical and 
quantum many body system, such as t-J model, Hub- 
bard model, and so on. It can help us to predict some 
new quantum PT. 

This work was supported by the National Natural Sci- 
ence Foundation of China. 
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